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We report on studies of the composition of ultra high energy cosmic rays with the Pierre Auger Observatory. 
The detection of longitudinal air shower profiles with the fluorescence detector is described and the measurement 
of the average shower maximum as a function of energy is presented. Furthermore, mass sensitive parameters 
that can be obtained from the observatory's surface detector data are discussed. 



1. Introduction 

The Pierre Auger Observatory is a facility to 
study the origin of ultra high energy cosmic 
rays by a combined measurement of their flux, 
anisotropy and mass composition. 

The southern site of the observatory is located 
near the town of Malargue in Argentina and its 
construction has been completed at the end of 
2008. Even during its assembly, data have been 
collected in a stable manner since January 2004. 
The lateral densities of charged air shower parti- 
cles at ground level are measured with the surface 
detector (SD), an array of 1600 water Cherenkov 
detectors distributed on an area of 3000 km 2 . 
In addition, the longitudinal air shower develop- 
ment can be observed with the fluorescence de- 
tector (FD), consisting of 24 wide-angle Schmidt 
telescopes that overlook the atmospheric volume 
above the surface detector. 

The mass composition of cosmic rays above 
10 17 eV is a key observable to discriminate be- 
tween different models put forward to explain the 
softening of the cosmic ray energy spectrum at en- 
ergies between 10 18 and 10 19 eV (the so called 'an- 
kle', see [2] for a recent review). Moreover, a mea- 
surement of the ultra high energy cosmic ray com- 
position is essential to understand the nature of 
the flux suppression observed above 5T0 19 eV [1]. 

In this article we will explain, how the cos- 
mic ray composition can be inferred from air 



shower measurements with the Pierre Auger Ob- 
servatory. An estimate of the primary composi- 
tion from the measurement of the longitudinal air 
shower profiles with the fluorescence detector will 
be given in the next section, followed by the de- 
scription of the mass sensitive observables of the 
surface detector. 

2. Fluorescence Detector Measurements 

Within the Heitler model [4] of cosmic ray in- 
duced particle cascades, the maximum number of 
electromagnetic particles of a shower with energy 
E and mass A is reached at a depth of 



(X laax )=a + X (]nE-(lnA)) 



(1) 



where Xq denotes the electromagnetic radiation 
length in air. The details of the hadronic inter- 
actions feeding the shower are subsumed into the 
parameter a = A In 2 — Xq In (3JVch-Ed.it) , where 
A denotes the hadronic interaction length, is 
the charged particle multiplicity and -Ecrit the en- 
ergy threshold, below which ionization becomes 
more important than radiative processes. The 
change of A max per logarithm of energy is called 
elongation rate [3]: 
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Thus, given a prediction of a and 
shower simulations, the primary mass composi- 
tion can be estimated from the average shower 
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Figure 1. Example of a measured longitu- 
dinal profile, reconstructed independently by 
two fluorescence detectors: Los Leones (black 
dots, E=31±2 EeV, X max =733±18 g/cm 2 ) and 
Los Morados (open squares, E=33±2 EeV, 
V max =751±16g/cm 2 ). 



maximum (A max ) and its relative change with en- 
ergy from the elongation rate. In the following we 
will describe the measurement of these two quan- 
tities with the Auger fluorescence detector. 

2.1. Event Reconstruction 

In the FD composition analysis only events 
recorded in hybrid mode [5] are considered, i.e. 
the shower development must have been mea- 
sured by the fluorescence detector and at least 
one coincident SD station is required to provide 
a ground level time that is needed for a precise 
determination of the shower direction. The longi- 
tudinal energy deposit profile is reconstructed [6] 
from the light recorded by the FD using the flu- 
orescence and Cherenkov yields from [7] and [8]. 
With the help of data from the observatory's at- 
mospheric monitoring devices [9] the attenuation 
of the light on its way from the shower to the 
detector is corrected for and the measured light 
emission heights are converted to atmospheric 
depth. Fig.[T]shows an example of a reconstructed 
longitudinal profile of an air shower, that trig- 



Figure 2. Sketch of the field of view of a fluo- 
rescence telescope and three longitudinal profiles. 
Showers with their maximum above FOVmax or 
below FOVmin are not selected. 

gered two of the four Auger fluorescence detec- 
tors. The data is fitted with a Gaisser-Hillas 
function [10] to extrapolate the profile to depths 
outside the field of view of the fluorescence tele- 
scopes. The integral over this function gives the 
calorimetric energy of the shower and the depth 
where it is maximal is at AT max . From showers 
that were observed simultaneously by more than 
one fluorescence detector, the difference in the re- 
constructed shower parameters can be used to in- 
fer the resolution of the X max reconstruction. It 
is about 20 g/cm 2 . 

2.2. Field of View Bias 

The energy and shower maximum can only be 
reliably measured, if V max is in the field of view 
of the telescopes (otherwise only the rising or 
falling edge of the profile is detected). As illus- 
trated in Fig. [21 this may lead to a biased event 
selection, if both, events with a small or large 
V max , are not included in the analysis because of 
their poor reconstruction. The current Auger flu- 
orescence detectors cover an elevation range from 
ill = 1-5° to SI2 = 30° and therefore the ob- 
servable heights for vertical tracks are between 
R tan f^i < h v < Rta,nH,2, where R denotes the 
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Figure 3. Field of view bias. Left and middle panel: Air shower simulations at E=10 19 5 eV. Right panel: 
Data at 10 18 eV < E <10 18 25 eV. 



distance of the shower core to the fluorescence 
detector. That is, the farther away from a fluo- 
rescence detector a track is detected, the smaller 
becomes the observable upper slant depth bound- 
ary X(FOV m in). Similarly the lower slant depth 
boundary X(FOV max ) becomes larger for near 
showers. 

The absolute magnitude of the field of view bias 
depends of course on the, a priori unknown, dis- 
tribution of the depth of shower maxima of cosmic 
ray showers. It is instructive to study the order of 
magnitude of the field of view bias for simulated 
showers with a known composition. The expected 
A max distribution for protons with an energy of 
10 19 5 eV is shown in Fig. 3(a) The upper view- 



able slant depth limits are indicated as arrows for 
different zenith angles. For vertical tracks it is 
around 880 g/cm 2 at the southern observatory, 
which is obviously not sufficient to detect the full 
X max distribution. A convolution of the trun- 
cated A max distributions with the arrival direc- 
tions of isotropic cosmic rays (diV/d cos 9 oc cos 0) 
leads to an effective bias for the measurement of 
the mean of the distribution. As can be seen in 
the upper panel of Fig. |3(b)| this bias can be as 
large as 12 g/cm 2 . Due to the unknown com- 
position and the theoretical uncertainties of air 
shower simulations it can not be corrected for a 
posteriori, since there is no unambiguous correc- 



tion given a measured (A max ). 
The strategy followed in this analysis is therefore 
to ignore showers that have a too small viewable 
slant depth range. The rejection of showers with 
A(FOV m in) <950 g/cm 2 almost completely re- 
moves the bias on (X max ) in the simulation, cf. 
lower panel of Fig. 3(b) The fiducial volume cuts 
that ensure an unbiased (X max ) can be derived 
directly from data as illustrated in Fig. 3(c) The 
dependence of (X max ) on the slant depth bound- 
aries is studied and events are only selected in the 
region where (X max ) is constant. 

2.3. Results 

The (X max ) results presented here were derived 
from hybrid data recorded between December 
2004 and April 2007. In order to ensure a good 
A max resolution, the following quality cuts were 
applied to the event sample: The reconstructed 
A"max should lie within the observed shower pro- 
file and the reduced x 2 of a fit with a Gaisser- 
Hillas function should not exceed 2.5. Moreover, 
insignificant shower maxima are rejected by re- 
quiring that the \ 2 °f a linear fit to the longitu- 
dinal profile exceeds the Gaisser-Hillas fit \ 2 by 
at least four. Finally, the estimated uncertainties 
of the shower maximum and total energy must 
be smaller than 40 g/cm 2 and 20%, respectively. 
Moreover, the fiducial volume cuts explained in 
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Figure 4. (A max ) as a function of energy [14]. 



the last section are applied. 

The systematic uncertainties of the atmospheric 
properties, event reconstruction and mass accep- 
tance add up to a total estimated systematic un- 
certainty of (X max ) that is < 15 g/cm 2 at low 
energies and < 11 g/cm 2 above 10 18 eV. 
In Fig. 2] the mean X max as a function of energy 
is shown along with predictions from air shower 
simulations [11]. As can be seen, the measure- 
ment favors a mixed composition at all energies. 

A simple linear fit, yields an elongation rate of 
54±2 (stat.) g/cm 2 /decade, but does not describe 
our data very well ( X 2 /Ndf= 24/13, P<3%). 
Allowing for a break in the elongation rate at 
an energy leads to a satisfactory fit with 
X 2 /Ndf= 9/11, P=63% and D m = 71 ± 5 
(stat.) g/cm 2 /decade below E h = io 18 - 35 eV and 
.Dio = 40 ±4 (stat.) g/cm 2 /decade above this en- 
ergy. This fit is indicated as a gray line in Fig. [4] 

Due to the uncertainties of hadronic interaction 
at highest energies, the interpretation of these 
elongation rates is, however, ambiguous. Using 
QGSJetII the data suggests a moderate lighten- 
ing of the primary cosmic ray composition at low 
energies and an almost constant composition at 
high energies, whereas the EPOS proton elonga- 
tion rate is clearly larger than the measured one 
at high energies, which would indicate a transi- 
tion from light to heavy elements. Theses ambi- 



Figure 5. Comparison to (X max ) from other ex- 
periments [12, 13]. 



guities will be partially resolved by the analysis 
of the A max fluctuations as an additional mass 
sensitive parameter. 

A comparison with measurements of other ex- 
periments is presented in Fig. [5] Taking into 
account the individual systematic uncertainties 
of each experiment there is a reasonable agree- 
ment of the (A max ) values. It is, however, worth- 
while noting that HiRes uses a different definition 
of A max (see [13]) that is about 10 g/cm 2 shal- 
lower than the one commonly used (depth of non- 
parametric shower maximum). If this is taken 
into account, the elongation rates of Auger and 
HiRes show a discrepancy in the X max -scale of 
about 30 g/cm 2 . 

3. Surface Detector Measurements 

As we have explained in the last section, the 
longitudinal air shower development is a good 
tracer for the primary composition that can be 
measured very precisely with fluorescence tele- 
scopes. However, these telescopes can only take 
data in moonless nights, whereas the surface de- 
tector has a duty cycle of almost 100% and cor- 
respondingly, the available event statistics are 
about a factor ten higher. Because of this, and 
the obviously independent systematic uncertain- 
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Figure 6. Simulated air shower signal in a water 
Cherenkov detector 

ties, it is worthwhile to study ground level ob- 
servables that are sensitive to the primary com- 
position as well. 

Since the amount of electromagnetic particles 
at ground depends for a given energy on the dis- 
tance to the shower maximum it depends also 
on the primary mass. Moreover, the number 
of ground level muons evolves differently with 
energy for different primary masses. From the 
Heitler model one expects 

oc E A 1 ' 13 (3) 

where p depends on the multiplicity and inelas- 
ticity of hadronic interaction (and is thus, as the 
above a, subject to theoretical uncertainties). 
Since the water Cherenkov detectors of the Pierre 
Auger Observatory cannot explicitly discrimi- 
nate between the electromagnetic and muonic air 
shower components on an event by event basis 
(but see Sec. 2]), a variety of 'indirect' experi- 
mental variables are currently investigated to re- 
late the surface detector data to the primary 
mass [22]: The signal rise time, its asymmetry 
and the signal shape analysifl 

1 For yet another composition sensitive parameter, the 
shower front curvature, see [21]. 



Figure 7. Correlation of rise time and X t 



3.1. Risetime 

For each event, the water Cherenkov detectors 
record the signal as a function of time (FADC 
traces) . This provides a possibility to distinguish 
the muonic and electromagnetic component, since 
the former travel in almost straight lines through 
the atmosphere, whereas the latter undergo mul- 
tiple scattering on their way to ground. There- 
fore, the two components have different path 
lengths and arrival times at ground. This is il- 
lustrated in Fig. [SJ where an example of a simu- 
lation [15, 16] of the surface detector response is 
shown. The rise time t\ji [17] is the time it takes 
until the cumulative signal grows from 0.1 to 0.5. 
It is a measure of the muon to electron ratio in 
a tank and depends on the primary mass, angle 
of incidence 6 and distance to the shower core 
r. The r and 9 dependence is marginalized by 
a comparison to the average rise time (£1/2(0, ?")) 
and the pull Aj = (t 1/2 - (£1/2(0, r)))/o-(£i/ 2 ) is 
averaged over the stations of an event to yield 
(Aj) [19]. This quantity could in principle be 
compared to air shower simulations. However, as 
previous studies showed, the number of muons 
measured in Auger can not be reproduced by 
these calculations [18]. Therefore, instead of a 
direct comparison of (A^) to air shower Simula- 
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Figure 8. Asymetry of station signals (0=54°, 
E=55±2 EeV) 



tions, it is related to the shower maximum. For 
this purpose, the subset of events that have been 
detected by both, the fluorescence and surface de- 
tector, are used to establish the relation between 
(Aj) and X max . As can be seen in Fig.[7J there is a 
linear correlation between the two variables, with 
which the surface detector data can be calibrated 
to yield X max ((Ai}), with which it is possible to 
measure the elongation rate with surface detector 
data. 

3.2. Risetime Asymmetry 

The SD stations measure one stage of the 
shower development for near vertical events. In 
case of inclined showers, however, considerably 
different shower ages can be observed, depend- 
ing on whether the station is up- or downstream 
of the incoming shower direction. This is illus- 
trated in Fig. [HI where the signal traces of two 
stations (both at the same shower plane distance 
r) are shown. The upstream trace is rather broad, 
which is compatible with a large fraction of elec- 
tromagnetic particles. The downstream trace is 
narrower indicating that most of the electromag- 
netic component has been attenuated and the sig- 
nal is dominated by muons. The corresponding 
asymmetry [20] in the rise time is thus a mea- 
sure of the longitudinal development of the ratio 



Figure 9. Distribution of jumps for events at a 
given energy, zenith angle and distance to the core 
range. 

of electromagnetic to muonic particles. The am- 
plitude of the asymmetry changes with the zenith 
angle, i.e. distance to the shower maximum. The 
angle at which it reaches its maximum can be 
used as an estimator for the primary composition. 

3.3. Signal Shape 

Since a muon deposits much more energy (typ- 
ically 240 MeV) in a water tank than an electron 
or photon (about 10 MeV), spikes are produced 
over the smoother electromagnetic background in 
the FADC time traces. Therefore, muons can be 
identified by searching for sudden variations, in 
the signal from one FADC bin to the next (cf. 
Fig [5]). The expected distributions of variations 
for purely electromagnetic and muonic traces can 
be fitted to the measured distribution as shown 
in Fig. [9] In that way, the muonic signal can be 
determined on a statistical basis with a resolu- 
tion of about 25% and the number of muons as 
a function of energy can be compared to predic- 
tions from air shower simulations to estimate the 
primary composition. 

4. Conclusions and Outlook 

In this article we presented the measurement of 
the elongation rate from data collected with the 
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fluorescence telescopes of the Pierre Auger Obser- 
vatory. When compared to predictions from air 
shower simulation the (A max ) data favors a mixed 
composition at all energies. This measurement 
will be soon updated with larger statistics and an 
analysis of the fluctuations of the shower maxi- 
mum. At lower energies the systematic uncertain- 
ties will be reduced by additional telescopes that 
cover the field of view from 30 to 60 degree [23]. 
A variety of surface detector observables are sen- 
sitive to the mass composition. These parame- 
ters have a somewhat worse mass resolution than 
I maJl which is, however, outweighed by the large 
event statistics collected by the ground array. 
Due to the hybrid design of the Pierre Auger 
Observatory, the correlation between the lateral 
and longitudinal shower parameters can be stud- 
ied and additional muon detectors will soon allow 
for an event-by-event measurement of the muon 
content of a shower [24]. 
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